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Abstract 

A thermodynamic stability criterion for the spontaneous breaking of translation 
invariance in many particle systems is derived. It simply requires the positive 
character of the wavevector dependent dielectric function, generalising the same 
condition for macroscopic dielectric constants. Its application to the Single Mode 
Approximation (SMA) for the description of collective modes of the v = 1/m 
Laughlin states is considered. The results indicate that the stability condition 
is violated by the SMA for all the relevant states m = 3, 5, 7, 9 in a wavevector 
neighborhood of the magneto-roton minima. These conclusions are in qualitative 
agreement with similar results obtained from the solution of the Bethe-Salpeter 
equation at v = 1/3 for both composite fermions and phenomenologically de- 
scribed electrons in the Laughlin state. 
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The stability properties of quantum many body ground states are of great interest 
in condensed matter theory. Important questions such as the occurrence of sponta- 
neous breaking of translation, and other, symmetries are closely linked with issues of 
stabilty. However, the complicated nature of the exact Green function equations for 
these problems usually makes extracting information about stability difficult. A partic- 
ular important example includes the determination of the conditions for the formation 
of the Wigner crystal at low electron gas density, both with and without a magnetic 
field §, §. 

In the present work, we discuss the derivation of a requirement to discuss the stabil- 
ity of a translationally invariant many particle system with respect to the spontaneous 
breaking of this symmetry. In short, this result is a generalisation of the usual posi- 
tive value condition on the dielectric constant H to the case of a general wavevector 
dependent dielectric function in a translationally invariant many body state. 

We expect this constraint may also be helpful in rigorously predicting the breaking 
of translation invariance which could be necessary to explain short range order in liquids 
that have fully invariant hamiltonians [[|. 

We apply the criterion to check the so called Single Mode Approximation (SMA) 
calculations of the collective mode excitations of the Laughlin ground states for 1/m 
values of the filling factor ||. 

The results indicate that the stability rule is not satisfied in a neighborhood of the 
magneto-roton minima for the states corresponding to m = 3, 5, 7, 9. The instability 
appears to be stronger at lower filling factors as is natural to expect. 

These conclusions are in qualitative agreement with our previous results || , [0 . In 
H, the collective mode dispersion for composite fermions at v = 1/3 was calculated 
numerically using the methods in ||. This work predicted an instability region in 
qualitative agreement with the one obtained here. In place of a magneto-roton mini- 
mum, the squared frequency crossed the zero frequency the momentum grew 
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from zero, and later became positive at higher momenta. 

In , we also considered the collective mode dispersion by using a phenomenolog- 
ical ansatz for the exact one electron Green function. It was also assumed that because 
the mean spatial separation between electrons is of the order of the expected size of 
the collective mode, screening effects in the interaction kernel of the Bethe-Salpeter 
equation were small. A region similar to the one found in |J and in the present work 
was obtained. 

It should be emphasized that the presence of an instability for these fractions does 
not necessarily mean that the Laughlin or composite fermion descriptions of these 
systems should be limited. These instabilities could be necessary for completing the 
picture given by the composite fermion description in which complementary flux quanta 
are associated with fermions in spatially wide approximate one particle Wannier states. 



This point of view was advanced in previous works , [ 1(| , ||TT| by one of the authors 
(A.C). However, its validity needs to be supported by the evaluation of ground state 
energies competitive with the ones corresponding to the Laughlin or composite fermion 
states. 

Moreover, other developments in the literature seem to support the spontaneous 
breaking of continous translation invariance in FQHE ground states. One which argues 
for a lower energy in periodically perturbed variants of the Laughlin wavefunction at 



v = 1/3 is discussed in | I2| . The presence of higher that expected peaks in numerical 
evaluations of the static structure factor was also treated in f[3" |. Another interesting 
result argued for the compatibility of the quantized Hall effect with the presence of 
crystalline structure in Quantum Hall crystals ||14|| . More recently [p~5f , the possibility 
of an instability in the composite fermion description has been suggested based on 
calculations of the effective mass for these particles. 



Special consideration requires the numerical evaluations for small systems ||16|| . It is 
clear that such calculations are essentially exact. It is also true that these evaluations 
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predict a ground state separated from the first excited level by a gap in qualitatively 
agreement with the value predicted by the SMA. Therefore, we conclude that numerical 
results for small systems do not support any strong spontaneous breaking of translation 
symmetry in the ground state. However, numerical calculations do not seem to be able 
to rule out the presence of a weak charge density or correlation periodic structure 
which could affect the gap at various filling factors through natural Bragg reflection 
processes. The highest particle number considered up to now is near twelve. With such 
a small number of particles, weak periodic properties present in the ground state in 
the thermodynamic limit, could be hidden by boundary problems. It is possible that 
the calculated energy gap between the ground state and the first excited level could 
also be viewed as arising from a weakly periodic state in the thermodynamic limit. 

From the experimental side, it is well known that finding traces of even the existence 
of the Wigner crystal phase with its strong charge density modulation is difficult. 
Therefore, the presence of a weak modulation in the properties of the ground state 
seems to be, at present, undetectable. 

In the next section, the stability criterion is given, and in Section 3 its application 
to the SMA is considered. 

1 Stability Condition 

Let us consider the following hamiltonian describing 2D-electrons with coulomb inter- 
actions 



H = / d 2 x^ + (x)0^(x) + / d 2 x(e^ + (x)^(x) - p )<j>{x) 




(1) 



The \1/ fields satisfy the usual equal time commutation relations 



+ 



(2) 
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which become in the first Landau level approximation of the QHE problem 



J + 



n (x-y), 



(3) 



for the specific fractional Hall effect system to be considered below . In (3), the n 
function is the projection operator onto the first Landau level. In this case, the free 



first Landau level restriction makes it a constant energy shift. A uniform background 
charge density p is assumed as usual to assure the stability of the system without 
breaking translation invariance. A macroscopic dielectric constant e is also considered 
for the medium in which the electronic system is imbedded. 

The 2D-electron gas interacts with an electrostatic field 4>(x) which is created by 
arbitrary but small fluctuations in the charge density of the background Pf(x). The 
field <f>(x) is therefore linked to Pf(x) through 



Assuming that the 2DEG is at a finite temperature determined by the parameter 
j3 = (fcT) -1 , the total internal energy of the electronic system plus the fluctuating 
electrostatic field is given by 



where the statistical mean is considered in terms of the Bloch density matrix p and the 
partition function Z 



particle kinetic energy term determined by O in (1) will be absent as usual because the 




(4) 




(5) 




p = exp(--), 
Z = Tr(p). 



(6) 



(7) 



(8) 



After using the expression for the entropy [pj]] 




< H > +kTlnZ, 



(9) 
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the total Free- Energy of the electronic gas plus electrostatic fluctuations can be written 
in the following form 

T[T,<f ) ] = U-TS = T e [T, 0] + J dx 3 ±- d4>{x)d<t>{x), (10) 

where the electronic contribution Y e is given by 

T e [T,0] = -±l n [Tr(eM-HP))\- (H) 

The starting point of our analysis will be the condition that the work needed to 
be done at constant temperature to produce the arbitrary background charge density 
fluctuation Pf(x) should always be positive for the system to be stable. Therefore, the 
second derivative of the Free-Energy functional over the field <fi should be a positive 
definite kernel. Taking the second derivative of (10) gives the following relation 

where 

6(x 3 ,y 3 ) = 5(x 3 )5(y 3 ). (13) 

The product of two 5 functions in the coordinates orthogonal to the 2DEG plane 
appears because the functional derivatives over are taken in 3D-space while the 
electron gas is confined to the x 3 = plane. 

To evaluate the electronic contribution in (12), it is convenient to use the Matsubara 
approach by introducing temperature Heisenberg field variables 



£ 4 ) = exp(x 4 H)ty(x) exp(-£ 4 i7), (14) 

ty+(x,x 4 ) = exp(x 4 H)^ + (x)exp(-x 4 H). (15) 

In terms of the particle density operator 

p(x, X4) = \l/ + (x, X4)ty(x, £4), (16) 
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the second derivative of the electronic contribution to the Free-Energy at vanishing 
values of the field (where the mean value of the density is homogeneous) can be 
expressed in the following form 

s$mm = -\ t dxt t dm K T ( ep(f ' Xi)M ^ Vi) - p2 ^ > 5(X3? m) 

= - dx A < ep(x, x 4 ).ep(y, 0) - p Q > S(x 3 , y 3 ) 
Jo 

rP . . 

= — dx4 < eexp(x4H)p(x).eexp(—x^H)p(x) > S(x 3 ,y 3 ) 
Jo 

+p 2 PS(x 3: y 3 ). (17) 

The energy eigenfunctions of the hamiltonian (1) satisfy 

H\n>=e n \n> (18) 

and a completeness relation for them can be introduced after the density operator p(x) 
in (17). This expression becomes 

<5 2 r e e 2 5(x 3 ,y 3 ) ^ exp[(-e m /3) - exp(-e n /3) 



54>{x)54>{y) Z n f^ 



-n <-m 



e 2 5(x 3 ,y 3 ) 



(3 < 0\p(x)\0 >< 0\p(y)\0 > +p 2 J8(x 3 ,y 3 ). (19) 



Z 

Let us consider now the zero temperature limit. As the temperature approaches zero 
the second term in (19) tends to cancel the constant contribution which is proportional 
to Pq since the partition function Z — > 1 (a non degenerate ground state is assumed) 
and the mean value densities are constant being balanced by the background jellium. 
In the remaining sum, only terms having one of the indices n or m zero will survive if 
the ground state energy is nondegenerate. Therefore, (19) reduces to 

S 2 T C , ,^ 2Re{< 0\p(x)\n X n\p(y>\0 >] 

After substituting (20) in (12), the second derivative of the total Free-Energy in 
the zero temperature limit takes the form 



5(x 3 , y 3 ) £ e 2^[<0\p(x)\n><n\p(y)\0>]^ ^ 

n^O €n £ 



We restrict in what follows the fluctuation charge densities to those that are non 
vanishing only in the 2DEG plane x 3 = 0. In this case the associated electrostatic field 
can be expressed in the Fourier representation 

00?) = / </>$) exp(ik.x-\k\\x 3 \), (22) 

— * — * 

where (fi(k) is any complex function of the momentum variables k in the plane satisfying 

m* = <j ) (-k). 

The above property follows from acting with the 3D-laplacian on (22) which pro- 
duces Poisson equation 

V 2 0O?) = - J J0p2\k\ (j>(k) exp(^.f) S(x 3 ) = -±np f (x), (23) 

where Pf(x) is the assumed arbitrary surface charge density fluctuation. 
Using this class of potentials allows us to write the quadratic form 

5 2 T 



5(j)(x)5(j)(y) 

; (2\k 
(2n)~ 



where the function K is given by 

y(t> V hM^ (25) 

and p£ is the 2D-Fourier transform of the density operator. 

Since the quadratic form (24) must be positive the following condition for the 
stability of the electronic gas arises 

■-WE^>». (26) 

eh\k\ v n ^o e n — e o 

This relation can be expressed in terms of the electron density p = N/V, the 
dynamical structure factor 

5 ft U>) = Jr E I < °\Pk\ n > - ( e « " e o)A), (W > 0) (27) 



and the interaction potential 

»(E) = 

ere 

resulting in the final expression 

l-2pv(k) / — S(k,w) >0. (28) 

JO w 

2 Analysis of the SMA 

In this section, the condition (28) will used to considered the Girvin, MacDonald and 
Platzman (GMP) single mode approximation (SMA) for collective modes in the FQHE 
ground states f|. For this purpose, it is convenient to introduce the dimensionless 
wavevector 

q = r Q k 

2 

and measure the energy in units of — where the magnetic radius is given in terms of 



the magnetic field intensity by r a = y ]|§|- 

As discussed in ||, the SMA is characterized by a dynamic structure factor of the 
form 

S(q,w) = s(q)5(w-A(q}), (29) 

where s is the static structure factor and A(q) is the collective mode energy as a 
function of the magneto-roton momentum. 
In terms of the polarizability parameter 

= |§, (30) 
and the electron density filling factor v: 

v 

P 



the stability condition for the approximate FQHE states described by the SMA can be 
written as 

y - Mo) > o- (si) 
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In order to investigate the validity of (31) for the SMA as applied to the main series 
of the v — 1/m Laughlin states, we used the data for the polarizability parameter 
obtained in || for m=3,5,7,9. The results indicate that the bound (31) is violated in 
a finite neighborhood of the maxima for the polarizability parameter as a function of 
the wavevector for all these states. The instability regions for v = 1/3 and v = 1/5 are 
shown in Figs. 1 and 2 as the zones of the curves lying above the straight lines. These 
lines are plots of the first term in (31) after multiplying by the factor m in each case. 
Similar results are valid for the v = 1/7 and v = 1/9 states. It can be seen that the 
instability is stronger for lower filling factors, as should be expected. 

In conclusion, a criterion that the static dielectric response must satisfy for an 
interacting many body theory to be translational invariant was presented. From it, 
indications of an instability in the static response of the Laughlin states within the SMA 
approximation were obtained. These results suggest the existence of states different 
than the SMA variational states which could break continuous translation invariance 
and have energies lower than the Laughlin wavefunctions. Such wavefunctions could be 



similar to those discussed in [12| which we suspect are related to the special Hartree- 
Fock states discussed in ||, jjlj and ||11|| . Alternatively, improvements of the original 
variational calculation of Girvin, MacDonald and Platzman may be possible. These 
problems will be considered in future work. 

The relation obtained here may also be useful for other problems where the stabil- 
ity of an homogeneous many body system under spontaneous breaking of translation 
symmetry may be an issue. 
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Figure Captions 

• Fig. 1 A plot of the SMA polarizability parameter as a function of the wavevec- 
tor for v — 1/3 Laughlin state. The momentum region in which the stability 
condition is not satisfied corresponds to the zone in which the curve is over the 
straight line. 

• Fig. 2 The plot of the SMA polarizability parameter for the v = 1/5 Laughlin 
state. The instability region is indicated by the wavevector neighborhood for 
which the curve passes over the straight line. 
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